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Abstract. This paper provides specific results on the Igusa local zeta function 
for the curves a;" In addition to specific results, we give an introduction 

to p-adic analysis and a discussion of various methods which have been used 
to compute these zeta functions. 



1. Introduction 

This paper provides specific results on the Igusa local zeta function for the curves 
f{x,y) = + J/™ where m and n are positive integers. We compute the following 
integral directly: 

Z{s)= f Ix"" +y"'\'dxdy 

for Re(s) > 0, where | • | represents the p-adic absolute value on the field of p-adic 
numbers, and dx denotes the Haar measure on normalized so that the measure of 
the p-adic integers Zp is 1. Igusa [5], [6] showed that this local zeta function has 
a meromorphic continuation to the whole complex plane and is, in fact, a rational 
function of i = p^*. We write Z{s) = Z{t). In addition to these specific results the 
paper gives an introduction to the more general methods and techniques of p-adic 
analysis and introduces Igusa's p-adic stationary phase formula and the resolution 
of singularities which are two other methods that have been used to compute local 
zeta functions. For more general computations of local zeta functions for forms of 
various types, please see [8], [M], [12], [4], [3]. 
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local zeta function. 
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2. P-ADIC NUMBERS. 



In this paper, we will let p stand for any fixed prime number, p e {2, 3, 5, 7, ...}. 
We can construct thep-adic numbers, Qp, from the rational numbers, Q, in exactly 
the same as the way in which we construct the real numbers, R. 

To do real analysis, we construct the real numbers from the rational numbers by 
defining the real numbers as the limits of all Cauchy sequences of rational numbers. 
In this way we insure that all Cauchy sequences of rational numbers converge to real 
numbers, or, in other words, that the real numbers form a complete field. Hence, 
this construction is called completing the rational numbers with respect to the usual 
absolute value, | ■ |oo = \/x^. When constructing Qp, we use a different absolute 
value, the p-adic absolute value. Thus, different sequences of rational numbers are 



Definition 1. The p-adic absolute value. A rational number, a, has a p-adic 
absolute value, \a\p, such that: 



The quantity ordpa is called the order of a and it is the highest power of p dividing 
a. In other words, if a = • ^, where a and /3 are integers that are not divisible 
by p, then n — ordpO. 



For example, if p = 5 then |75|5 = 15^ • 3|5 = S^^, |34|5 = |5° • 34|5 = 1, and 



Definition 2. A p-adic Cauchy sequence. A sequence of rational numbers, {a„}, 
is said to be a p-adic Cauchy sequence if given some e > 0, there exists an integer 
N such that for i,j > N, \ai — aj\p < e. 

Two p-adic Cauchy sequences, {a„} and are said to be equivalent if ja; — 

b'ilp — >■ as i — > oo. For example, if p = 5 then {a„} = {1, 1, 1, ...} and {bn} = 
{1 + 5, 1 + 5^, 1 + S'', ...} are equivalent since |1 — (1 + 5*)|5 — > as i oo. 

Definition 3. The field of p-adic numbers, Qp. The p-adic numbers, Qp, are 
defined to be all equivalence classes of p- adic Cauchy sequences. 

Every equivalence class of Cauchy sequences in the real numbers has a standard 
representation, which we know as the decimal expansion. Likewise, it can be shown 
that every equivalence class of Cauchy sequences in the p-adic numbers has a unique 
representative of a certain form, which is called its p-adic expansion. The form of 
the decimal expansion does not determine a unique real number; for example, the 
sequences {1,1,...} and {.9, .99, .999, ...} are equivalent Cauchy sequences which 
both represent the real number 1. For a detailed explanation of the p-adic numbers 
and their expansions, see [TT] . 

Theorem 1. |llj The ring of p-adic integers, Zp. Every equivalence class, a, of se- 
quences of rational integers has exactly one representative p-adic Cauchy sequence, 
{an}, where a„ € Z, of the following form: 



Cauchy. 





{a„} = {ao, flo-l-aip, oq aip -\- a2p'^ , ...} 

where ai — oq + a\p + a2P^ + ... + a^p* and aj £ {0, 1, ...,p — 1}, for < j < i. 
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These equivalence classes with their unique representative form the p-adic in- 
tegers, denoted Zp. It is convenient to think of such an equivalence class as 
its representative expansion and we denote a G Zp as an infinite series: a = 
ao + aip + a2P^ -f .... This series representation of the p-adic integers is analogous 
to the decimal expansion for real numbers, and we often write a = 00.010203.... 
We define the norm of a € Zp, to be if a represents and in the non-zero case 
\oi\p = p~°''''p("), where ordp{a) is the exponent of the lowest power of p with non- 
zero coefhcient in the p-adic expansion of a. Any rational integer, n, is a p-adic 
integer because it is the limit of the constant sequence {n, n, ...} and has a termi- 
nating p-adic expansion (to signify the repeating sequence), for example, for p = 5, 
405 = l-5-fl-52 + 3-53 = 5(1 + 1 -S-l-S-S^) and |405|5 = 5"^ Thus, the ordp{a) 
is the highest power of p which divides the p-adic expansion of a. However, there 
are more p-adic integers then the rational integers, for example, if p = 5 then 1/6 
is also an integer because it can be written asl4-4-5 + 4-5^ + ....(To check that 
this expansion represents 1/6 multiply it by 6 = 1 -I- 1 • 5 to get 1.) However, not 
every rational number is a p-adic integer. The number 36/5 = 1/5(1 -I- 2 • 5 -I- 1 • 5^) 
is not a 5-adic integer; it is an element in Q5 of absolute value 5. It can be shown 
that Qp, as we defined it above, is equivalently defined as the field of fractions of 
Zp. Hence, we can think of elements in Qp as: 

a ^ ao + QlP+Q2P^ + ... ^ (ord(a)-ord(/3))y 

P feo + 61P + 62p2 + ... 
where a € Zp, /? G Zp — {0} and m is a unit in Zp. The units in Zp are the 
invertible elements and they have p-adic absolute value 1. In other words, u is a 
unit if the first digit in its p-adic expansion is not 0. For example, we can find the 
5-adic expansion of the number 30/225 by long division (where we think of as 
5-|-4-5-|-4-5^ + ... and borrow from the right if we need to) as follows: 

30 5(1 -H- 5) 1 2 ^3 4 

225 = 5^(4 + 1. 5) = 5(^ + '-' + ^'' +^'' + 

In the representation above, 6/9 = 2/3 is a 5-adic unit and we can sum the conver- 
gent geometric series above to show that we have the correct expansion. 

°° 1 c2(n+l) 2 

4+l-5+3-52+l-53+3-5*+... = 4+5(1+3-5) V52' = 4+5(16) lim 5— 77 

n— »oo 1 — 5 3 

i=Q 

For more examples of numerical computations with the p-adics see [T] , [TT] , [T3] . 

The p-adic integers, Zp, form a ring with base p addition and multiplication, 
where carrying is done to the right. The set of p-adic integers with first digit 
equal to are divisible by p, hence denoted pZp, and form the maximal ideal in 
Zp. The ideal p^Zp is the set of p-adic integers with first two digits equal to 
and p'^Zp is the set of p-adic integers with first e digits equal to 0. Hence, we see 
that Zp D pZp D p^Zp D ... D p'^Zp D ... D {0} is a descending chain of ideals. 
The units in Zp are denoted by Zp — pZp since they are elements in Zp but not 
in pZp. The units can be written as a disjoint union of the p — 1 cosets of pZp: 
Zp — pZp — UaGfi 2 + P'^p)- The p-adic integers themselves are an infinite 

disjoint union: 

00 00 
Zp = [] p'^ Zp - p'^+i Zp = ]J p^ (Zp - pZp) 



The p-adic numbers can be visualized as an infinite system of concentric circles 
representing the ideals p^Zp for all integer values of e and corresponding to the 
inclusions {0} C ... C p^Zp C pZp C Zp C p~^Zp C .... Since is contained in all 
the ideals, it lies at the center of all the circles. If the order of a p-adic number a is 
e then a falls into the p^Zp circle but not into the 'p^^^'L^ circle and hence wc think 
of it as an element in the annulusp'^(Zp— pZp). In this way, we see that the disjoint 
union above for the f)-adic integers is really a union of annuli, and, more generally, 
Qp, itself, is a union of infinitely many disjoint annuli, each annulus containing the 
elements in the field of the same absolute value: 

oo oo 

Qp= ]J p«Zp-p«+iZp= ]J p^(Zp-pZp). 

e=— oo e= — oo 



3. The definition of the Igusa local zeta function and p-ADic 

ANALYSIS. 

The Igusa local zeta function, Z{s)^ associated to a polynomial in the ring 

of polynomials with n variables and coefficients in Z, j{x) = f{xi,X2,---,Xn) € 
Z[xi,X2, --jXri], is defined for a complex variable s, where Re{s) > 0, as : 

Z{s) = I \fix)\;dx. 

As |/(a;)|p is always a power of p~^, we can think of Z{s) as a function oit = 
and write Z{t). 

Before we continue, we must explain what we mean by integration in this setting. 

Definition 4. Measure of the set E. Let be a union of sets Ui, U2, ■■■Um of the 
form Ui = tti +p^'Zp (analogous to intervals in real analysis), where Oj G Zp, Ni 
is any non-negative integer, and 1 < i < m. We define the measure of the set E, 
■m{E), to be the integral over E (written m{E) = J^dx) such that the following 
four properties hold: 

1. m(S) = Jj^dx>Q and m(0) = 

2. Translation invariance. If a G Zp, m{a + E) = J^_^^ dx = J^dx = m{E) 

3. If UifMJj = 0, then m{Ui\JUj) = J^.^^. dx = J^, dx+ J^, dx = m{Ui)+m{Uj) 

4. m{Zp) = J^^dx = l. 

This measure is called the Haar measure on Qp normalized so that m(Zp) = 1 
and it is known to be the unique measure to satisfy these four properties. 

Proposition 2. The measure of the ideal p^Zp is p~^ for any integer e. 

Proof. We break up Zp into its p*^ cosets modulo p'^Zp and use properties 2,3, and 
4 of the measure. Each coset is of the form bo + bip + feaP^ + ... + b'^~^p'^~^ + p'^Zp 
and, as there are p choices for each bi, we get p'' cosets. Hence, 

1 = / dx = / dx = p'^ j dx. 
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From this argument, we see that m{p'^Zp) = p □ 

Remark 1 . This proposition gives a change of variables formula which we write as 
d{p^x) = \p'^\pdx = p~'^dx. In fact, by the same proof the following, more general 
proposition is true: If £' is a measurable set of the type given above then the 

measure of the set p'^E is p~'^m[E). 



4. The Igusa local zeta function for the polynomial f{x) — 



Break Zip into a disjoint union (as above) of annuli where l^-'^lp is constant. 
Thus, 

Z{t)^ I \x''\;dx = Y, I \x''\ldx. 

Now, making the change of variables x = p'^y for y £ Zp — pZp and dx = p~'^dy. 
We see that 

m -flp-'f \{p'y)Xdy = f^pi-^-^^)^ f dy= y— ^ / 

e=0 "'Zp-pZp JZp-pZp ^ P JZp- 

As the units are a disjoint union of their p — 1 cosets modulo pZp, we have that 

/ dy^Y.1 dy={p-l)j dy = I ~ p-\ 

And we have that the Igusa local zeta function for f{x) — x^ is precisely 

= i_p-i-Ns = 1 -p-HN- 
This integral is extremely well-known and appears in Tate's Thesis |15) . 

5. Local zeta functions and their Poincare series. 

As defined above, the Igusa local zeta function associated to the polynomial f{x) 
in n variables with coefficients in the integers ( or equivalently in the p-adic integers) 
has a Poincare series which contains arithmetic information about the cardinality 
of solutions of f{x) in the ideals p'^Zp. 

Before we consider the Poincare series, let's look more closely at Z{t). The zeta 
function is really an infinite sum of the measures of all points in that are mapped 
into each annulus of Zp by the polynomial /(x). Hence, 

m = / \I{x)\ldx = ^m({x e T;\f{x)=p-^u}) 

•^^P e=0 

where u € Zp—pZp. From this expression for Z{t), it is clear that Z{1) = 1 and that 
Z(0) = m{{x\f{x) = u\). If we know this first measure and f{x) is a homogeneous 
polynomial whose partial derivatives vanish modulo p only at 0, then by Hensel's 
Lemma we can compute the measures of the other sets and find the rational function 
for Z{t) in terms of this first number. The problems in computing the Igusa local 
zeta function come when the polynomial has complicated singularities. 



The Poincare series which is associated to the Igusa local zeta function is a 
generating function for the number of solutions of /(x) modulo p'^ for all non- 
negative integer e. It is defined as: 

oo 

P{t)=J2\Ne\p-"'t^ 
e=0 

where n is the number of variables in f{x), \Ne\ is the cardinality of the set: = 
{a e Zp — p'^Zp |/(a) = mod p'^}, and A^o = 1- Borewicz and Shafarevic conjecture 
the rationality of P{t) in [5], p. 47. Igusa proves this conjecture in [5], [B] by showing 
that Z{t), and hence P{t), is rational. 

Proposition 3. [5], [6] The local zeta function associated to a polynomial is related 
to the polynomial's Poincare series as: 

Z{t)^P{t)-\{P{t)-l). 



Proof. If / ^(p'^Zp) is the pull-back (or inverse image) of p'^Zp then 

/-i(p^Zp) = {xeZ;|/(x) = Omodp^}= [] (a + p'^Zp). 

Hence, m(/~^(p"^Zp)) — m(/^^(j3'^+^Zp)) is the measure of the set of points in Zp 
that map to the annulus p'^Zp — p^+^Zp. With this in mind, we see the following: 



Z{t) = / \f{^)\ldx 



= y\ / p'^'dx 

^ J/-i(p=Zp)-/-i(p-+iZp) 

OO 

= Y^p-<^^[m{f-\p<^Z,)) m{f-\p<^+'Z,))]. 

e=0 

We remember that /"^(p'^Zp)) = UaGW {o.'^P'^'^p)- The number |iVe| is the number 
of centers a e Zp — p'^Zp such that /(a) = Omodp*^. The measure of a disk 
a + p^Z^; is since /„+z„ da: = Ja,+z, ^xi J^^+^J^^^- J.^+zJ'^n = P""". 

Thus, m(/-i(p'=Zp)) = I^Veb"^''". Now we see that 



OO 



Z{t) = J2 \Ne\ip-''tr - J2 We+llp-^^+'^'^t^ 
e=0 e=0 

= P{t)^t-\P{t)-l). 

□ 



6. Direct computation of the Igusa local zeta for + y™. 

As an example to familiarize the reader with p-adic analysis, we will calculate 
the Igusa local zeta function for f{x,y) ~ + y™ directly. We will show that 

Z{t) = ^ [ p^LWnJ-Ji^J + p-Un/mi-j^nj^ 



^ [{p - If - (|iv(0)| - i)]p-Hi - p-H) + (|jv(0)| - i)ip^ - p)p-H 

{l-p-H){l-p~^tT^) 

First we break Zp up into its infinitely many annuli lJ^of'^(^p ~ P^p) and write 
it as 



J'^p e=0 •^P°(^P~P^p) 

Similarly, we break up the region of integration for the other variable, and make 
the change of variables x p^x and y — >■ p^y, to get that 



oo oo „ 

Z{t) = Y.Y.P~^P~^ / If*"'^" + P^'^y^'lpdxdy. 

/=0e=0 i(Zp-pZp)2 

Note that x and y are now both units in (Zp —pZp). To compute the integral above, 
we break it into three partial integrals. We consider the three cases: 1) ne > mf, 

2) mf > ne, and 3) m f = ne. 

In the first case where ne > mf, we have a factor of p"^-^ in the integrand. In 
this case, we can write the partial integral as 

f>Oe>mf/n J CZp-p'LpY 

Since y is a unit, the integral is equal to the measure of (Zp — pZpY which is 
{1 — p~^Y and we have that 

f>Oe>mf/n 

Now, 

oo Ymf/n\ 

(1) E p-' = T.p~'- E p-" 

e>mf/n e=0 e=0 

We need to sum from e = to the greatest integer in mf /n which we will 
denote by \mf/n\. This sum depends upon the congruence class of / modulo 
where (to, n) is the greatest common factor of to and n. If f = , . where 



( 

fc is a positive integer then we sum from e = to , . . If f = , "'^ , + 1 then 
™//" = J^+^f^ we sum from e = to + [m/n\ . If / = then 
m/Zn = T^^fc +2m/n and we sum from e = to + \ 2m/n\. If / = + 7' 
where < j < ^ " n — 1 then mf In = -nr^ + jm/n and we sum from e = to 

Therefore, 

E ^'-^= E E p-''= E ^ (i-p-i) — - 

e=0 j=0 e=0 j=0 ^ ^ ' 



Using formula ([1} above, we have that the inner sum in Ii (t) is 



e>mf/n j—0 

Using this sum, Ii (t) becomes 



fe>0 j=0 

Similarly in the second case, when mf > ne, the second partial integral becomes: 
Therefore, we see that 

Finally in the third and most difficult case where mf = ne, 

hit) = V p-™//"-/r^' / \f{x,y)\'dxdy 
f>o. ^ez Az.-pz,)= 



E(P"^^'^)' / \f{x,yWdxdy 



fc=0 
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(Zp-J9Zp)2 



(|iV(0)|-l)(l-p-i)p-2t 



(1 - p (™"^") i (".") ) "'(Zp-pZp) 

and we will show that 

m= f \fix,yWdxdy=[ip-l)^-i\Nm-l)]p , 

where f{x,y) = a;™ + y", and p does not divide both m and n, |iV(0)| is the 
cardinality of the set N{0) = {(x, y) € Fp x Fp | + y" = modp} = A^i + 1, 
and Fp is the finite field with p elements. 

First proof using counting argument 

(Zp -pZp)2 = UZoii^^y) G U'\\fi^,y)\' = t"}, where i - p-^ So, 

OO 

m^J2*' m{{x,y)eU' I \f{x,y)\^^t^}. 

k=0 



Hence, 



I'sit) = lim V t*^ m{{x, y) € (c^ + cip + ... + Cr-ip"" ^) +p''Zp 



r— ^oo 

fc=0 



x{dQ + d\p + ... + dr-ip''~^) +p"Lp 

I Co, do e F;, a, di e Fp for i > 0, \f{x, y)\' = t''}. 

This interpretation means that 

r 

I^{t) = lim card{{x, y) e (co + Cip + ... + Cr-ip^'^) x (do + dip + ... + dr-ip^'^) 

/c=0 

Now since = if and only if ordp{x"^ + y") = fc. If wc let \ord = 

k\r denote the cardinality of the set {{x,y) £ (cq + Cip + ... + Cr-ip^^^) x {do + 
dip+ ...dr-ip'^~^) I |/(a;,y)'' = i'^}. Then we see that I^{t) = limr^oo X]fe=o l^''^ ~ 
k\rt''p-'^'-. 

What is \ord = A;|r? 

\ord = 0\r= card{{x,y) e (cq + cip + ... + c^-ip''"^) 

X {do+dip+ ... + dr-ip''~^) I ordp(x'" + y") = 0}. 

Now ordp{x"^ + y") = if and only if p does not divide (co + Cip + ...Cr-ip'^~^)"^ + 
{do + dip+ ... + dr-ip'^~^)'^ which happens if and only if p does not divide Co^ + do" 
by the binomial theorem. Thus, ci,...,Cr-i and di, ...,dr-i do not influence the 
condition ordp{x'^ + y") =0 and are each free to take any value in Fp. With this 
in mind, we see that 

lord = 0\r = card{{co,do) G (f;)^ | co^+do") = modplp^C-D = [(p_i)2_(^(o)_i)]p2r-2_ 

\ord = l\r = {ix,y) € {cQ+cip+...+Cr-if^^)x{do+dip+...+dr-ip''^^) \ ordpix'^+y"^) = 1} 

Now, ordp{x"^ +2/") = 1 if and only if p divides x"^ + y" and does not divide 
_j_ yTi -y^rjiipji happens if and only if p does divide c™ + dp but does not divide 
(co + cip+ ... + Cr-ip^^^)™ + (do + dip+ ... + This last condition means 

that p^ does not divide (cq + cip)™ + (do + dip)" which by the binomial theorem 
implies that does not divide co™ + mcoCip + do" + ndodip. Now, if p divides 
Co™ + do", then co™ + do" = p"^, for some e > 1 and we have that ordp{x"^ +y'^) = 1 
if. p^ does not divide p"^ + mcoCi|5 + ndodip. This last condition implies that p does 
not divide p''^^ + mcoC\ + ndodi. Now, Co and do are units, so if to is a unit, then 
p^/x™ + y" implies that ci = (— p*^"^ — ndodi)TO~"'^CQ ^ mod p. Thus, di can be any 
element of Fp, and ci is determined by di. So, there are p choices for (ci, di) G Fp 
such that p^ divides a;™ + y", and therefore (p^ — p) choices such that p^ does not 
divide + y". Similarly, if n is a unit, then there are again (p^ — p) possibilities 
for (ci,di) such that p^ does not divide x™ + y". Finally, 

|ord=l|r = card{(co,do)e (F;) I a;™ + y" = Omodp} 

X card{{ci,di) € (F^) | a:'" + y" ^ modp^} 

X card{(c2,...,Cr_i,d2,...,dr_i) e (Fp)^^'-^)} 

= (7V(0)-l)(p2-p)p2'-*. 
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In general, for A: > 1, \ord = k\r = card{{co,do) G (F*)^ + = mod 
p} ■ card{{ci,di) e | a;™ + y" = mod p^} . . . card{{ck-i,dk-i) |a;™ + = 
modp*^} • card{{ck,dk) \x"^ + y" ^ modp''+^} hence 

lord = k\r = (NiO) - - 

Now, we can compute our integral. 



I^{t) = \im p-^''y\ord = k\rt'' 

k=0 

r 

= lim p-^^liip - if - {N{Q) - l)^''-^ + y((A^(0) - 

r 

= ((p - 1)2 - (Ar(0) - + lim y ((A^(0) - l)p^{p' - p)p'^''-^-''h^+^) 

oo 

= ((p- 1)2 - (iV(0) - l))p-2 + (iV(0) - l)(p2 -p)p-4f^p-'=t'= 



fe=0 



^2 r.rra. .^^„-2 , - 1)(1 - P'^-^t 

1-p-H 



{{p-ir-{N{o)-i))p-'+ 



In the case where m = 2 and n = 3, -/V(O) = p, so 



7^(t) = ((p-l)2_(p_l))p-2 + 



2, (p-l)(l-p-l)p-2f _ (p-l)(p-2)p-2+f(l-p-l)p- 



1 -p-ii 



1-p-H 

Second proof using a change of variables 

We will compute I^{t) = J^^ ^ |x™ + y"!" dxdy using another method. Write 
the integral as a sum over its cosets modulo p and count those points for which 

_|_ yU -g g^j^j jg jjQ^ g jnodulo p. 

+y2|« dxdy 



a;™ + y"|* dxdy = 



(Zp-pZp)2 



(a,6)eF;xF;-^(«>'')+pZ^ 

1)2 _ (1^(0)1-1)) 



dxdy 



+ i\N{0)\-l)p- 



Ka+px)™ + (6 + py)"|' dxdy 



/Z2 

P 

Expand (a + px)™ + (& + py)" to get |(a™ + ma^-^a; + ■ • ■ + p'"x™) + (6" + 
nb"~^py + . . . + p"y")| = |p(c + ma'^~^x + nb'^~^y + p{. . .)| where a™ + 6" = cp 
for some c € F*. There is a measure-preserving transformation (x,y) — >■ (x*,y*) 
such that X* = yi(x,y) = ma™~^x + nb"'~^y + p{. ■ .) and y* = g2{x,y) = y. The 
Jacobian of this transformation is 



det 



dgi 


9gi 




dx 
dg2 


dy 
992 


= det 


dx 


dy 





dgi 

dx 



1 

^m— 1 



(a, b) = ma'" 



If ma™"^ = mod p, then let gi{x, y) = x and 52 = ma"'~^x + nb"~^y + p{. . .) so 
that the Jacobian becomes n6"~^ ^ mod p since p does not divide both m and 
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n. For a general description of this measure-preserving map see [8,9]. Now, we see 
that \f{a+px,b + py)\^ = \p{x* + c)\^ where dxdy = dx*dy*. Returning to the 
integral above, 

m = {{p-if-{\Nm-i))p-^+{\Nm-i)p-H( \x*+c\' dx*dy* 

(|7V(0)|-l)p-2t(l-p-i) 



((p-l)^-(|Ar(0)|-l))p- 



l-p-H 



Hence, Z{t) = Ii{t) + l2{t) + /3(t) and we have computed the local zeta function 
in terms of [7V(0)|. Below we compute |-/V(0)| for our curve. 

Theorem 4. For any prime p the number of solutions to a;™ + = in Fp, which 
we have called |A^(0)| above, is 

1 + {p — l)gcd{m.,n,p — 1) if p = 2 or 

iford2{p— 1) > min{ord2{m),ord2{n)} 



\N{0)\ = 



otherwise. 



Proof. The case p = 2 is trivial. 

If ord2{p — 1) > min{ord2{m),ord2{n)} then w.l.o.g. assume ord2{p — 1) > 
ord2{m). The set Bt = {a:*|a; G F*} of t powers forms a subgroup of order 
{p — l)/gcd{p — l,t) in F*. Let Z = Bm,{^Bn then Z is a subgroup of order 

9cd{\B^\,\B^\) = Since ord2{p - I) > ord2{m) = Y, 2^+^ 

divides p — I and there is an clement x in F* which generates the cyclic subgroup 
of order 2^+^ in F*. Since x generates this subgroup, x^^ = —1 and as is 
odd we see that x'^ = x^^ ^ = (-1)^ = -1. 

Allowing for the point (0, 0), we have that |A/'(0)| = 1 + |{(o, 6) e F* x F*|o™ + 

6" = mod By the existence of the particular clement x above, we have that 
|A/-(0)| = 1 + |{(a,fo) e F; X F;|o™ = 6" moduli since for each a where a™ = c 
there is a unique ax where [ax)™' = —c. As a™ = 6" implies that a™ and 6" are in 
the subgroup Z above, we can rewrite \ 'N(fS)\ = 1 + Yl,z<^z \{^'^^ ^) € F* x F*\a"' = 
= z mod p}| = 1 + E^ez{a ^ F;|a" = z)\ x \{b € F;|6" = z}\. For each z 
in Z, \{a e F;|a™ = z}\ = gcd{p - l,m) and \{b € F;|6" = z}\ = gcd{p - l,n). 
Therefore, |A''(0)| = l + \Z\ x gcd{p-l,m)gcd{p-l,n) = l + {p-l)gcd{p-l,m,,n). 

If ord2{p — 1) < min{ord2{m), ord2{n)} and c = 2°^'^^^p~^\ we know that c 
divides m and n. Now if we have (a, b) ^ (0, 0) such that a™ + 6" = mod p then 
= o" + 6" = a™ + 6™ implies that (a")'= = -(6")'=. As 5 ^ 0, we have the 
element y = a"(fe^^)'" such that y'^ = —1. Thus y^"^ = 1 and the order of y must 
divide 2c = 2°'''^2(p-i)+i jf ^]^g order of y < 2c then the order of y would divide 
c implying the impossible y° = 1. Therefore, the order of y must be equal to 2c 
but as 2c does not divide p — 1, we have a contradiction and the only solution to 
a™ + &" s mod p is (0,0). Hence, in this case \N{0)\ = 1. □ 

7. Computation of the Igusa local zeta function for f{x, y) = a;" + 

USING THE p-ADIC STATIONARY PHASE FORMULA (SPF) 
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Theorem 5. A p-adic stationary phase formula [10 . The local zeta function of a 
polynomial f{x) in n variables with coefficients in Tip is: 

Zit) = (p" - |A^(0)|K" + (|7V(0)| - \S\)p-^H^-^ + ( \f{x)Ydx 

where \N(f))\ is the cardinality of the set of x in for which f{x) = mod p, \S\ 
is the cardinality of the set of singular x in N{0) such that all partial derivatives of 
f are modulo p at these points, and S is the set of all x in which are congruent 
to vectors in S mod p. 

Corollary 6. To compute Z{t) — dxdy, SPF must be applied 

-r^ + - 1 times. 



Proof. For the curve x™' + y^\ SPF will terminate when the original integral appears 
as the singular integral. We may assume that m < n. For if to = ti, SPF will 
terminate after one application. The first application of SPF results in 

/ Ix'^ + y"'\'dxdy = h + h +P^^ f Ip'^x'^ +p"'x"'\' dxdy 

Jzl 

Where Ii and I2 are known in terms of |iV(0)|. Taking out the factor of p™ the sin- 
gular integral becomes p^^t™ ^2 |p"~'"a;" + y™!'' dxdy. Applying SPF to this 
integral produces the known terms and a new singular integral p~^ /„2 + 
pmym^B (j^j.(iy_ 'J'q proceed, we must determine whether n — to > to or n — m < to. 
In other words, we have to know which is bigger n or 2m. 

The string of SPF applications ends when instead of > or < we get = in which 
case, all powers of p factor out of the integrand and the singular integral becomes 
/z xz 1^" + y"^\^ dxdy. Taking a multiple of Z{t) to the other side, we are now 
able to solve for the rational function. 

With each application of SPF we get a new inequality of the form kn < Im or 
kn > Im. Suppose we may assume this to be the inequality for the SPF of 

+p"y™|" dxdy, 

Zp X Zp 

SO since kn ^ {I ~ l)m < to, the integral reduces to 

And our next application of SPF produces the integral 

|p"j;" + p'^-'^"?/™!" dxdy 

Zp X Zp 

which can be simplified after determining which is greater, (fc + l)n or Im. As we 
can see, each application of SPF increases either fc or Z by 1, and since SPF must 
be applied exactly once more when k = and I = , " x and fc = ^ = 1 at the 

start, the total number of applications of SPF will be , " s + / " ^ — 1. □ 
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8. Computation of the Igusa local zeta for x" + y"* using resolution 

OF singularities. 



Resolve f{{x, y j) ~ x"^ + to find the zeta function 
Z{t)= [ \f{{x,y))\' dxdy = J2 [ \foh\'\h*{dxAdy)\ 

JZpXZp j-)^j-)^Jh-^{ZpXZp) 

Where h* is the pull back of h and h~^{Zp x Zp) = Dj for some i. Let /i : M 
such that 



and 



Therefore, 



So, 



Mte,r/i)) = (2/,^). 



j^-'i-Uy. ")") — / ('''I'^iyi) condition described below 
U ^y)) ^ condition described below 



Resolve f{{x, y)) = xl{xi + yf). 



h{{x2,y2)) = (yi,— ) 
2/1 

h{{i2,m)) = {xi, — ) 

Xl 

Therefore, 

{X2y2,X2) 



h ^{{xi,yi) 
So, 



(6,6%) 



f((x v)) - <; ^1^2(2:2 + 2/2) 



Now, we resolve f{{x,y)) = X2y|(a;2 + 2/2) 



Therefore, 



h{{x3,y3)) = {X2,^) 
Mte,%)) = (2/2,-) 

2/2 



((^2,2/2))- I ^6^3(1 + ^3) 



Neither of these curves need to be resolved further. Now we must find the D^. Since 
X and y are in Zp, |a;|p < 1, \y\p < 1, so either |||p < 1 or |||p < 1. Therefore 

\xi\ < l,|2/i| < 1 and |6| < Ij*?!! < 1- This means that (xi.yi) e Z^ and 
{£,i,r]i) G Zp X pZp. Therefore, Di = Zp x pZp and the same argument holds for 
D2, D3, while Di = Zl. 
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Next we calculate the change in measure due to the change of variables in Di. 
dx Ady = d{^i,rji) A d^i 

= 6 drji A d^i 

= -^1 d^i A dr]i 

Now, calculate the first piece of the zeta function over Di . 



^ 1 - p-2i2 

Next we calculate the change in measure due to the change of variables in D2. 

dx A dy = dxi A d{xiyi) 
= Xi dxi A dyi 

= 6 ^6 Ad(^2'?2) 

= (.2 d^2 A dr]2 

Now, calculate the second piece of the zeta function over D2. 



Next we calculate the change in measure due to the change of variables in D3. 
dx Ady = X2y2 d{x2y2) A dx2 





= P 



-1 I-P-' 
1 - p-3t^ 



= -xly2 dx2 A dy2 

= -(6%)^6 d{^3Vs) A d^3 

= ^hl d^3 A drja 



Now, calculate the third piece of the zeta function over D^. 
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Next we calculate the change in measure due to the change of variables in D4. 



dx f\dy = -xlx^uz dxs A 1^(0:3^3) 

= -x^ys dx3 A dys 
Now, calculate the fourth piece of the zeta function over D4. 

Ziit) = [ \4yU^ + ys)\' ^3^2/31 dxsdys 

V 

= f \x,\^'+Uxs f |2/3|'^+'|2/3 + ir rfy3 



[1 

[2: 

[4 

[6: 
[7: 

[10: 

[11 
[12: 

[13: 

[14 



The Igusa local zeta function, Z{t), is the sum of Zi,Z2,Z3,Z4, or 

(l-p-H)(l-p-5t6) 
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